Calorimetric measurement of quantum work 
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To define the work performed on a driven quantum system in a physically sound way has turned 
out to be a truly non-trivial task, except in some special cases of limited applicability pHllj. This 
topic has been in a focus of intense research recently in the attempts to generalize the classical 
fluctuation relations (T34I31 into the quantum regime. Here we propose and demonstrate that 
a calorimetric measurement gives both a theoretical and experimental tool to test the Jarzynski 
equality (JE) and other fluctuation relations in a quantum system, and to analyze the distribution 
of dissipation in them, based on the very principle of conservation of energy. We focus on an 
experimentally feasible two-level system, a superconducting Cooper pair box [T^ - f20( ] subject to 
Landau-Zener interband transitions [2l| - [2a |. Because of the small heat capacity and weak relaxation 
to the phonon bath, the calorimetric measurement on the electron gas (resistor) turns out to be a 
very feasible experimental method. 



Consider the basic setting for non-equilibrium fluctua- 
tion relations [l2, 14 1: initially the system is coupled to a 
bath and in thermal equilibrium with it. Thereafter it is 
driven by a control field in a non-equilibrium fashion. We 
are interested in the distribution of work performed on 
the system in repeated experiments with identical driving 
protocol among the realizations. In order to secure ther- 
mal equilibrium in the next realization, the interval be- 
tween the repetitions is ideally infinitely long. We mainly 
consider the protocols where the drive signal returns to 
its initial position or to an equivalent one with equal free- 
energy at the end. Our main statement in this respect 
is then that, due to the energy conservation, the work 
W done on the system under these conditions is equal 
to the heat Q dissipated to the environment in the time 
interval from the beginning of the driving till the end of 
the equilibration period. With this premise we can avoid 
the formidable task of analyzing the work itself during 
the quantum evolution, and characterized by difficult-to- 
analyze time-orderings of the operators involved. The 
calorimetric principle can, at least in principle, be ex- 
tended to more general protocols by bringing the system 
back to the initial position adiabatically at the end of the 
protocol. 

Before going into the detailed analysis of the calori- 
metric method, let us first define the quantum system 
and present results in certain limiting cases for the work 
distribution, which can be obtained by known methods. 
We focus on a quantum two-level system, shown in Fig. Q] 
and characterized by an avoided level crossing under the 
operation of the control parameter q. For simplicity, we 
consider a symmetric set-up, where the maximum energy 
spacing between the ground (g) and the excited (e) level 
is -Emax in the beginning and the end of the drive, and 
Emin at the minimum in the middle. We normalize q such 



that it obtains the values =Fl/2 at the beginning and the 
end of the drive, respectively, and assumes a value q = 
at the avoided crossing. Suppose for the basic illustration 
that the system evolves in a unitary way from the initial 
state: this represents very weak coupling to the environ- 
ment. In other words, the system is initially thermalized 
by coupling it to a bath with inverse temperature j3, and 
decoupled from the bath for the duration of the driven 
evolution. This corresponds to a common situation for 
considering the fluctuation relations 13, 12|. The value 
of work in each realization is then determined by the 
internal energy stored in the system during the opera- 
tion. There are two possibilities: the system makes a 
Landau-Zener (LZ) transition around q = between the 
two states with probability Plzj or it does not make the 
transition at probability 1 — Plz- If the system makes 
the transition, work ±E max is done on the system (+ for 
g — > e transition, and — for the e — > g transition) by the 
end of the driving, otherwise the performed work van- 
ishes. The distribution of work W can then be written 
as 

P (W) = (1 - p° gg )P LZ S(W + E max ) 

+ (1 - Plz) S(W) + P ° gg P LZ 6(W - E max ), (1) 

in agreement with the results in Refs. 0, 0]- Here 
6(E) refers to the Dirac delta function, and p gg = 



(1 



c ) 1 is the initial thermal population in the 



ground state. We may also state that the same distribu- 
tion as for work applies for heat in the sense discussed 
above, since after an infinitely long time, for even weak 
coupling, the system has relaxed back to the thermal 
equilibrium state by exchanging energy with the environ- 
ment. Examples of work distributions (with broadening) 
are shown in Fig. [T] (b) at a few bath temperatures. As 
for the quantum systems evolving unitarily from a ther- 
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FIG. 1: A two-level system with avoided crossing between 
the ground (g) and the excited (e) state. In (a) the control 
parameter q is swept from left to right and the system is ini- 
tially in the ground state and evolves either in the ground 
state (probability 1 — Plz) or makes a Landau-Zener transi- 
tion (probability Plz)- (b) Examples of the distribution of 
work in repeated sweeps under the conditions that the system 
starts in thermal equilibrium at various temperatures, and it 
is not coupled to the environment during the evolution. The 
distributions correspond to equation {TJ with artificial broad- 
ening of the peaks, and for Plz = 0.3, and /3-B max = 0.1,0.3 
and 1.0 (red, blue and black curves, respectively). The peak 
at W = — -Emax arises from e to g LZ transitions, the peak at 
W = is for adiabatic evolution, and that for W = +i?max 
arises from g to e transitions, as indicated by the diagrams 
within the (b) panel. 



mal state, the JE holds naturally for the two level system. 
Indeed we see immediately that 



(e-f> w ) = 



dw P (wy 



1 



(2) 



for the distribution of equation ([1}. Here the brackets 
(•} refer to averaging over the experimental realizations. 
Equations ([T]) and ([2]) simply confirm the earlier under- 
standing of JE in a closed quantum system. The true 
challenge, which can be addressed by the here-proposed 
calorimctric detection, is, however, the measurement of 
work in an open quantum system coupled to the environ- 
ment also during the driving period. 



Although our proposal and arguments are quite gen- 
eral, we specify next to a particular physical system, a 
Cooper- pair box (CPB), where concrete results about the 
distribution of work, and a description and analysis of the 
calorimetric measurement are immediately possible and 

~ 2C 



can be implemented experimentally. The CPB [IE 
consists of a superconducting island connected to a su- 
perconducting lead by a Josephson tunnel junction. The 
system is described by the circuit scheme in the inset of 
Fig. [5] and it is characterized by a voltage source V g , 
coupling gate capacitance C g , a Josephson junction with 
energy Ej and capacitance Cj. We denote Cs = C g +Cj. 
Resistor R forms the dissipative environment of the box. 
In the regime e = Ej/(2E C ) < 1, where E c = 2e 2 /C s is 
the charging energy of the box, we can treat the CPB as 
a two-level quantum system [20j. Denoting with |0) and 
|1) the states with zero and one excess Cooper-pairs on 
the island, respectively, the Hamiltonian of the system 
reads 



H 



-2fc?(|l><l|-|0)<0|) 



■(|1)(0| + |0}(1|), (3) 



2 Thus the system 
Ej at q = 0, see 



where q = C g V g /(2e) is the normalized gate voltage. 
We assume a linear gate ramp q(t) = —1/2 + t/r over 
a period r = 1/q starting at t = 0. The energy gap 
separating the ground and the excited state of the sys- 
tem is given by AE = 2Eq \J q 2 + e 
passes a minimum energy gap E, 
Fig. Q](a), where it makes a transition with the probabil- 
ity Pl z = exp(— ™ ^ CT ) according to the standard LZ 
model [2ll423| . The value of the energy gap at the ends 
of the trajectory is E max = E C \/l + (Ej/E c ) 2 ^ E C - 

We proceed further by exploiting the calorimetric prin- 
ciple to an open driven two-level system. The detailed 
dynamics of the CPB under the driven evolution with LZ 
transitions in the dissipative regime is beyond the focus 
in the present work. Instead, we consider the two level 
system in the incoherent regime, where energies of the 
order of the minimum gap Ej or smaller are irrelevant. 
Then the system dynamics is governed by a master equa- 
tion with rates between charge states |0) and |1) given 
by 



T ± (AE) = -EjV(±AE), 



(4) 



where AE = 2Ec<q is the energy dissipated (absorbed) 
in the + (— ) tunneling event and exchanged as a photon 
between the resistor and the system, and 



V(E) 



1 

2?r^ 



cftcxp(J(t) + -Et) 



(5) 



is the corresponding probability density of the system to 
emit energy E to the environment 24 1. The standard 
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correlation function in equation ([5]) is given by 



4 f c 

7T Jo 



OJ Rq 



{ coth(-/?fru;)[cos(cjt) — 1] — isin(wf)}. (6) 

Here, i?Q = ?i/e 2 and the relevant impedance of the CPB 
reads 



9?eZ t (w) = 



R 



err 



1 + (i? cff C*erfw) 5 



(7) 



where i? eff = A 2 i?, C cff = (A -1 - 1)C S , and A = C S /C E 
is the coupling parameter. To obtain the distribution of 
work, our task then is to keep track on the net energy 
deposited to the resistor in each ramp. 

We have analyzed the distribution of work numerically 
by stochastic simulations tracing the paths, assuming 
rates of equation (|4]). Results of such simulations are 
shown in Fig. [2] For all these distributions JE is numer- 
ically valid within 0.1 %, and it should naturally be valid 
identically since V(E), and thus the ± rates in equa- 
tion (0| obey detailed balance (25|. These simulations 
correspond to the pumping trajectory described above, 
with parameters given in the Fig. [2] caption. The dis- 
tributions have several features that can be explained by 
simple arguments as follows. First of all, the peaks in 
(a) and the steps in (b) at W = ±Ec, and at cor- 
respond to the peaks arising from excitation/relaxation, 
and from the adiabatic dynamics of the unitary evolu- 
tion in equation ([T]), respectively, and they get gradually 
weaker towards faster relaxation (i.e., towards increasing 
R). In general, the mean of the work is positive, i.e. the 
evolution is dissipative. Furthermore, we have checked 
numerically that the Crooks equality [14| is valid in the 
form p(—W) = e~" w p(W) for these distributions. 

Next we discuss the calorimctric measurement of the 
energy deposited on the resistor. The electronic system 
of the resistor forms the environment that is then cou- 
pled to a super-bath formed, e.g. by phonons to which 
these electrons couple by standard electron-phonon cou- 
pling. This leads to the following sequence once an en- 
ergy quantum AE is emitted/absorbed by the two-level 
system due to =p transitions, see Fig [3] We assume the 
standard hot-electron regime l27j . (i) The temper- 
ature T of the resistor (electrons) changes abruptly by 
AT = ±AE/C at the moment of the photon exchange; 
here C is the heat capacity of the electron gas in the resis- 
tor, (ii) The electronic temperature of the resistor starts 
to relax back towards the equilibrium To = l/(fcs/?) at 
the rate governed by the thermal conductance G that 
is determined, e.g. by electron-phonon coupling, (iii) 
New energy pulses may occur during the evolution. In 
this evolution the deposited energy (heat) to the bath 
can be obtained generally from Q = CT + GAT, where 
AT = T— Tq. Note that the pulses can have either + or — 




FIG. 2: Numerical results of the distribution of work in a 
CPB. The circuit is shown in the inset of (b). In (a) we 
present the actual distribution p(W). The parameters of 
the system are Ej/k B = 0.1 K, E c /k B = 1 K, A = 0.3, 
bath temperature fcs//3 = 0.2 K, r = 1 • 10~ 9 s, and 
R/Rq = 1/3, 1, 10/3, 10. Each distribution is based on 9- 10 7 
repetitions. In (b) we show the corresponding cumulative dis- 
tributions U(W) = f_ x p(W')dW with the same colouring 
of the lines as in (a). 



sign, i.e., the resistor may cool as well. In the considered 
example, Q(t) = JT AEi8{t — ti), where AEi is the mag- 
nitude of heat released in the ith absorption/emission 
event, and U are the time instants of these events. Note 
that unlike in the analysis leading to equations ([1]) and 
@, we allow here different values for AEi for transi- 
tions occuring along the driven trajectory. Integrating 
the equation for Q for the situation where the system is 
in thermal equilibrium in the beginning and at the end, 
i.e. T(0) = T(oo) = To, and assuming small temperature 
rise in the pulses, we have for the heat in each realization 



Q = G AT{t)dt. 



(8) 



To make full use of equation ©, the following practical 
issues need to be taken care of [27[ ■ The thermometer 
is first calibrated against the bath temperature under 
the equilibrium conditions (without applying the driving 
field), or one uses a primary thermometer that serves 
without calibration. The only remaining calibration is 
then the determination of G: this can be obtained by 
applying a known constant power P to the calorimeter, 
and measuring the temperature rise ST of the electrons 
in the thermometer, since G = P/ST. 

Before turning to the detailed practical realization and 
numbers, one further consideration is in order. Namely, 
biasing the gate capacitor in the ramp that we consider 
leads to a current that dissipates Joule power AEd in the 
resistor R. This dissipation can be estimated by elemen- 
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FIG. 3: Schematic presentation of the calorimetric measure- 
ment. The resistor in the middle emits (absorbs) heat to 
(from) the two-level system on the left upon excitation (re- 
laxation) shown by blue (red). The corresponding heat pulses, 
i.e. temperature of the electrons in the resistor vs. time, are 
shown on the right with the same colour conventions. We 
assume that the electrons are in the so-called hot-electron 
regime, where the internal (electron-electron) relaxation is 
faster than other time scales in the problem. 



tary circuit analysis with the result AEd = e 2 R/r. For 
typical values of r and R, for instance those in Fig. [5J 
this energy is at most of order O.OlEc- Since Ec is the 
scale of work performed, the Joule heating produces only 
a small error in the calorimetric measurement, even if no 
precautions are taken. 

Finally we analyze the feasibility of measuring calori- 
metrically the distribution of quantum work in a su- 
perconducting two-level system considered. The most 
straightforward realization would involve a small litho- 
graphic metal resistor, ~ (0.1 fim) 3 , on the chip [28j . 
The electronic heat capacity at To = 100 mK tempera- 
ture is C ~ 10~ 20 JK _1 , yielding a jump in temperature 
of Ec /C ~ 3 mK upon absorbing the relaxation heat at 
the end of the ramp, assuming a realistic Ec/ks = 3 
K. This is a sizable change of temperature (3%) which 
can be detected over a relatively long time of order 10 -4 
s that is determined by the equilibration of the resis- 
tor back to its initial temperature via electron-phonon 
coupling. The large sensitivity of the calorimeter is nat- 
urally due to the standard T n /Tp ~ 10 ~ 6 suppression 
of the heat capacity of the electron gas with Fermi tem- 
perature Tp ~ 10 5 K. If needed, longer relaxation times 
can be obtained by, for instance, suspending the resistor 
on a silicon nitride membrane [29(, which would make 
the detection even easier. A tunnel probe, in form of 
a superconducting lead connected to the resistor, i.e. a 
SIN junction, measured by RF reflection or transmission 
techniques [30( , is fully compatible with the requirements 
of the presented calorimetric detection of heat. 



[1] G. N. Bochkov, G. N. & Kuzovlev, Yu. E. General theory 
of thermal fluctuations in nonlinear systems. Sov. Phys. 
JETP 45, 125-129 (1977). 

[2] Kurchan, J. A Quantum Fluctuation Theorem, e-print 
cond-mat/0007360. 

[3] Tasaki, H. Jarzynski Relations for Quantum Systems and 



Some Applications, e-print cond-mat/0009244. 

[4] Mukamel, S. Quantum Extension of the Jarzynski Re- 
lation: Analogy with Stochastic Dephasing. Phys. Rev. 
Lett. 90, 170604 (2003). 

[5] Yukawa, S. A Quantum Analogue of the Jarzynski Equal- 
ity. J. Phys. Soc. Jpn. 69, 2367-2370 (2000). 

[6] Chernyak, V. & Mukamel, S. Effect of Quantum Collapse 
on the Distribution of Work in Driven Single Molecules. 
Phys. Rev. Lett. 93, 048302 (2004). 

[7] Allahverdyan, A. E. & Nieuwenhuizen, Th. M. Fluctu- 
ations of work from quantum subensembles: The case 
against quantum work-fluctuation theorems. Phys. Rev. 
E 71, 066102 (2005). 

[8] Engel, A. & Nolte, E. Jarzynski equation for a simple 
quantum system: Comparing two definitions of work. 
EPL 79, 10003 (2007). 

[9] Talkner, P., Lutz, E. & Hanggi, P. Fluctuation theorems: 
Work is not an observable. Phys. Rev. E 75, 050102(R) 

(2007) . 

[10] Campisi, M., Hanggi, P. & Talkner, P. Quantum fluctua- 
tion relations: Foundations and applications. Rev. Mod. 
Phys. 83, 771-791 (2011); 83, 1653 (2011) [E]. 

[11] Solinas, P., Averin, D. V. & Pekola, J. P. Work 
and its fluctuations in a driven quantum system. 
arXiv. 1206.5699. 

[12] Jarzynski, C. Nonequilibrium Equality for Free Energy 
Differences. Phys. Rev. Lett. 78, 2690-2693 (1997). 

[13] Bochkov, G. N. & Kuzovlev, Yu. E. Nonlincr fluctuation- 
dissipation relations and stochastic models in nonequi- 
librium thermodynamics. I. Generalized fluctuation- 
dissipation theorem. Physica A 106, 443-479 (1981). 

[14] Crooks, G. E. Entropy production fluctuation theorem 
and the nonequilibrium work relation for free energy dif- 
ferences. Phys. Rev. E 60, 2721-2726 (1999). 

[15] Crooks, G. E. Path-ensemble averages in systems driven 
far from equilibrium. Phys. Rev. E 61, 2361-2366 (2000). 

[16] Jarzynski, C. Nonequilibrium work relations: founda- 
tions and applications. Eur. Phys. J. B 64, 331-340 

(2008) . 

[17] Averin D. V. & Pekola, J. P. Statistics of the dissi- 
pated energy in driven single-electron transitions. EPL 
96, 67004 (2011). 

[18] Averin, D. V., Zorin, A. B. & Likharev, K.K. Bloch os- 
cillations in small Josephson junctions. Sov. Phys. JETP 
61, 407-413 (1985). 

[19] Bouchiat, V., Vion, D., Joyez, P., Esteve D. & Devoret, 
M. H. Quantum coherence with a single Cooper pair. 
Physica Scripta T76, 165-170 (1998). 

[20] Nakamura, Y., Pashkin, Yu. A. & Tsai, J. S. Coher- 
ent control of macroscopic quantum states in a single- 
Cooper-pair box. Nature 398, 786-788 (1999). 

[21] Landau, L. Zur Theorie der Energieubertragung. II. 
Phys. Z. Sowjet. 2, 46-51 (1932). 

[22] Zener, C. Non-adiabatic Crossing of Energy Levels. II. 
Proc. R. Soc. A 137, 696-702 (1932). 

[23] Shevchenko, S., Ashhab, S. & Nori, F Landau-Zener- 
Stiickelberg interferometry. Phys. Rep. 492, 1-30 (2010). 

[24] G. L. Ingold and Yu. V. Nazarov, in Single Charge Tun- 
neling, edited by H. Grabert and M. H. Devoret, NATO 
Advanced Study Institutes, Ser. B, Vol. 294 (Plenum 
Press, New York, 1992), p. 21. 

[25] Crooks, G. E. Nonequilibrium Measurements of Free En- 
ergy Differences for Microscopically Reversible Marko- 
vian Systems. J. Stat. Phys. 90, 1481-1487 (1998). 



5 



[26] Richards, P. L. Bolometers for infrared and millimeter 
waves. J. Appl. Phys. 76, 1-24 (1994). 

[27] Giazotto, F., Heikkila, T. T., Luukanen, A., Savin, A. M. 
& Pekola, J. P. Opportunities for mesoscopics in ther- 
mometry and refrigeration: Physics and applications. 
Rev. Mod. Phys. 78, 217-274 (2006). 

[28] Timofeev, A. V., Helle, M., Meschke, M., Mottonen, M. 
& Pekola, J. P. Electronic Refrigeration at the Quantum 
Limit. Phys. Rev. Lett. 102, 200801 (2009). 

[29] Vercruyssen, N., Barends, R., Klapwijk, T. M., Muho- 
nen, J. T., Meschke, M. & Pekola, J. P. Substrate- 
dependent quasiparticle recombination time in supercon- 



ducting resonators. Appl. Phys. Lett. 99, 062509 (2011). 
[30] Schmidt, D. R., Yung, C. S. & Cleland, A. N. Nanoscale 
radio-frequency thermometry. Appl. Phys. Lett. 83, 
1002-1004 (2003). 

Acknowledgements We thank J. Ankerhold, S. Gas- 
parinctti, D. Golubcv, F. Hckking, M. Mottonen, and 
J. Santos for discussions. The work has been supported 
partially by LTQ (project no. 250280) CoE grant and the 
European Community's Seventh Framework Programme 
under Grant Agreement No. 238345 (GEOMDISS). 



